Gravitational waves generated by second order effects during inflation 
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The generation of gravitational waves during inflation due to the non-linear coupling of scalar and 
tensor modes is discussed. Two methods describing gravitational wave perturbations are used and 
compared: a covariant and local approach, as well as a metric-based analysis based on the Bardeen 
formalism. An application to slow-roll inflation is also described. 



I. INTRODUCTION 

The generation of gravitational waves (GW) is a gen- 
eral prediction of an early inflationary phase [l[. Their 
amplitude is related to the energy scale of inflation and 
they are potentially detectable via observations of B- 
mode polarization in the cosmic microwave background 
(CMB) if the energy scale of inflation is larger than 
- 3 x 10 15 GeV SSSSli- Such a detection would be 
of primary importance to test inflationary models. 

Among the generic predictions of one-field inflation [tJ 
are the existence of (adiabatic) scalar and tensor per- 
turbations of quantum origin with an almost scale in- 
variant power spectrum and Gaussian statistics. Even 
if non-linear effects in the evolution of perturbations are 
expected, a simple calculation [8] , confirmed by more de- 
tailed analysis Q , shows that it is not possible to produce 
large non-Gaussianity within single field inflation as long 
as the slow-roll conditions are preserved throughout the 
inflationary stage. Deviations from Gaussianity can be 
larger in, e.g., multi-field inflation scenarios [1, [13] and 
are thus expected to give details on the inflationary era. 

As far as scalar modes are concerned, the deviation 
from Gaussianity has been parameterized by a (scale- 
dependent) parameter, /nl- Various constraints have 
been set on this par ameter, mainly from CMB analy- 
sis QJ] (see Ref. [TJ for a review on both theoretical 
and observational issues). Deviation from Gaussianity in 
the CMB can arise from primordial non-Gaussianity, i.e. 
generated during inflation, post-inflation dynamics or ra- 
diation transfer [14]. It is important to understand them 
all in order to track down the origin of non-Gaussianity, 
if detected. 

Among the other signatures of non-linear dynamics is 
the fact that the Scalar- Vector and Tensor (SVT) modes 
of the perturbations are no longer decoupled. This im- 
plies in particular that scalar modes can generate gravity 
waves. Also, vector modes, that are usually washed out 
by the evolution, can be generated. In particular, second- 
order scalar perturbations in the post-inflation era will 
also contribute to B-mode polarization (l5j or to multi- 



pole coupling in the CMB [16] . and it is thus important 
to understand this coupling in detail. 

In this article, we focus on the gravitational waves 
generated from scalar modes via second order dynam- 
ics. Second-order perturbation th eory has been investi- 
gated in various works [H EE EI Hfl HI HI H H HI 
and a fully gauge-invariant approach to the problem was 
recently given in Ref. [25| . Second-order perturbations 
during inflation have also been considered in Refs. d Hi| , 
providing the prediction of the bispectrum of perturba- 
tions from inflation. 

Two main formalisms have been developed to study 
perturbations, and hence second order effects: the 1 + 3 
covariant formalism [2^ ] in which exact gauge-invariant 
variables describing the physics of interest are first identi- 
fied and exact equations describing their time and space 
evolution are then derived and approximated with re- 
spect to the symmetry of the background to obtain re- 
sults at the desired order, and the coordinate based ap- 
proach of Bardeen [28| in which gauge-invariants are 
identified by combining the metric and matter pertur- 
bations and then equations are found for them at the 
appropriate order of the calculation. In this article we 
carry out a detailed comparison of the two approaches up 
to second order, highlighting the advantages and disad- 
vantages of each method, thus extending earlier work on 
the linear theory [2i|. Our paper also extends the work 
of Ref. [22[], in which the relation between the two for- 
malisms on super-Hubble scales is investigated. In partic- 
ular, we show that the degree of success of one formalism 
over the other depends on the problem being addressed. 
This is the first time a complete and transparent match- 
ing of tensor perturbations in the two formalisms at first 
and second order is presented. We also show, using an 
analytical argument, that the power-spectrum of gravi- 
tational waves from second-order effects is much smaller 
than the first order on super-Hubble scales. This is in 
contrast to the fact that during the radiation era the 
generation of GW from primordial density fluctuations 
can be large enough to be detected in principle, though 
this requires the inflationary background of GW to be 
sufficiently small [23| . 
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This paper is organized as follows. We begin by re- 
viewing scalar field dynamics in Section II within the 1 
+ 3 covariant approach. In Section III, we formulate 
the problem within the covariant approach followed by a 
reformulation in the coordinate approach in Section IV. 
A detailed comparison of the two formalisms is then pre- 
sented in Section V. In Section VI, we study gravitational 
waves that are generated during the slow-roll period of 
inflation. In particular, we introduce a generalization of 
the /jvl parameter to take into account gravity waves 
and we compute the three point correlator involving one 
graviton and two scalars. Among all three point func- 
tions involving scalar and tensor modes, this correlator 
and the one involving three scalars are the dominant Q . 
Finally, we conclude in Section VII. 

II. SCALAR FIELD DYNAMICS 

Let us consider a minimally coupled scalar field with 
Lagrangian density 1 

Ct = -V=g[lv a cj>v a <t> + v(<f>)] , (2) 

where V{<j>) is a general (effective) potential expressing 
the self interaction of the scalar field. The equation of 
motion for the field <f> following from is the Klein - 
Gordon equation 

V a V a 0-F'(0) =G, (3) 

where the prime indicates a derivative with respect to <f>. 
The energy - momentum tensor of 4> is of the form 

T ab = V a 0V 6 - g ab [\V c <t>V c 4> + V((f>)] ; (4) 

provided <j) ta ^ 0, equation ([3]) follows from the conser- 
vation equation 

V b T ab = . (5) 

We shall now assume that in the open region U of space- 
time that we consider, the momentum density V a is 
timelike: 

V Q 0V a </> < . (6) 

This requirement implies two features: first, <fi is not con- 
stant in U, and so {<j> = const.} specifies well-defined 



1 We use conventions of Ref. [3Qj. Units in which h= c = kg = 1 
are used throughout this article, Latin indices a, b,c... run from 
to 3, whereas Latin indices i,j,k... run from 1 to 3. The symbol 
V represents the usual covariant derivative and d corresponds 
to partial differentiation. Finally the Hilbert-Einstein action in 
presence of matter is defined by 

A = j dx^y^ij^R + C^,} . (1) 



surfaces in spacetime. When this is not true (i.e., <fi is 
constant in U), then by ((4|), 

V a </> = & T ab = -g ab V(<t>) => V = const., (7) 

in U, [the last being necessarily true due to the conser- 
vation law |5j))] and we have an effective cosmological 
constant in U rather than a dynamical scalar field. 



A. Kinematical quantities 

Our aim is to give a formal description of the scalar 
field in terms of fluid quantities; therefore, we assign a 4- 
velocity vector u a to the scalar field itself. This will allow 
us to define the dot derivative, i.e. the proper time deriva- 
tive along the flow lines: f a - b c ... d = u e V e T a ' b c ... d - 
Now, given the assumption ([6]), we can choose the 4- 
velocity field u a as the unique timelike vector with unit 
magnitude (u a u a = —1) parallel to the normals of the 
hypersurfaces {(f) = const.} [3l| 2 , 

u a = -i^W, (8) 

where we have defined the field ip = (f> = {— Va^V^) 1 / 2 
to denote the magnitude of the momentum density (sim- 
ply momentum from now on). The choice © defines 
u a as the unique timelike eigenvector of the energy - 
momentum tensor (j4|. 3 

The kinematical quantities associated with the "flow 
vector" u a can be obtained by a standard method [H,[3^. 
We can define a projection tensor into the tangent 3- 
spaces orthogonal to the flow vector: 

h ab = g ab + u a u b => h a b h b c = h a c , h ab u b = ; (9) 

with this we decompose the tensor V b u a as 

V b u a = V b u a - u a u b , V b u a = ±<dh ab + a ab , (10) 

where V a is the spatially totally projected covariant 
derivative operator orthogonal to u a (e.g., V a / = 



2 In the case of more than one scalar field, this choice can still be 
made for each scalar field 4-velocity, but not for the 4-vclocity 
of the total fluid. A number of frame choices exist for the 4- 
velocity of the total fluid, the most common being the energy 
frame, where the total energy flux vanishes (see 1331 for a detailed 
description of this case. 

3 The quantity if) will be positive or negative depending on the ini- 
tial conditions and the potential V; in general <f> could oscillate 
and change sign even in an expanding phase, and the determi- 
nation of u a by (0 will be ill - defined on those surfaces where 
Va0 = => tp = (including the surfaces of maximum expan- 
sion in an oscillating Universe). This will not cause us a problem 
however, as we assume the solution is differentiable and (J6j holds 
almost everywhere, so determination of u a almost everywhere by 
this equation will extend (by continuity) to determination of u a 
everywhere in U. 
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h a b ^bf', see the Appendix of Ref. [H[ for details), u a 
is the acceleration (u b u b = 0), and a ab the shear (a a a — 
c a bU h = 0). Then the expansion, shear and acceleration 
are given in terms of the scalar field by 



°ab = -^ _1 /l<a C /lb> d V c [V d 0], 

a a = -i/j^ 1 V a ip = -V' _1 (V a V' + u a ip) , 



(11) 
(12) 
(13) 



where the last equality in Eq. (fTI) follows on using the 
Klein - Gordon equation ([3]) . We can see from Eq. flB"]) 
that ip is an acceleration potential for the fluid flow [361 ] . 
Note also that the vorticity vanishes: 



■Jab = -h a c h b d \7 [d (V'~ 1 V c]( /.) = 



(14) 



an obvious result with the choice ([5]), so that V a is the 
covariant derivative operator in the 3-spaces orthogonal 
to u a , i.e. in the surfaces {<f> = const.}. As usual, it is 
useful to introduce a scale factor a (which has dimensions 
of length) along each flow -line by 



(15) 



where H is the usual Hubble parameter if the Universe 
is homogeneous and isotropic. Finally, it is important to 
stress that 



V a <£ = 



(16) 



which follows from our choice of u a via equation (|8|). a re- 
sult that will be important for the choice of gauge invari- 
ant (GI) variables and for the perturbations equations. 



B. Fluid description of a scalar field 

It follows from our choice of the four velocity ([5]) that 
we can represent a minimally coupled scalar field as a 
perfect fluid; the energy - momentum tensor (H|) takes the 
usual form for perfect fluids 



T ab = fiu a u b + ph ab 



(17) 



where the energy density (i and pressure p of the scalar 
field "fluid" are given by 



H = \^ 2 + V(<f>) , 

P = - v{4>) • 



(18) 
(19) 



If the scalar field is not minimally coupled this simple 
representation is no longer valid, but it is still possible to 
have an imperfect fluid form for the energy - momentum 
tensor [3l[ . 

Using the perfect fluid energy - momentum tensor (JTTJ) 
in ([5]) one obtains the energy and momentum conserva- 
tion equations 



lp 2 U a + ^aV = 



(20) 
(21) 



If we now substitute n and p from Eqs. (fT8| and (fi"9| into 
Eq. (|2H)) we obtain the 1+3 form of the Klein- Gordon 
equation J3J: 



4> + Q<j> + V'(<p) = 



(22) 



an exact ordinary differential equation for <p in any space - 
time with the choice ([5]) for the four - velocity. With the 
same substitution, Eq. ([2Tj) becomes an identity for the 
acceleration potential i\>. It is convenient to relate p and 
H by the index 7 defined by 



P = (7 - 1)/-^ & 7 



p + M ip 2 



li 



(23) 



This index would be constant in the case of a simple one- 
component fluid, but in general will vary with time in the 
case of a scalar field: 

7 V 
1=6(7-2-2—. (24) 

7 ip 

Finally, it is standard to define a speed of sound as 

2 - P 1 7 fOK\ 

c ^ = 7 - 1 -e 7 -' (25) 



C. Background equations 

The previous equations assume nothing on the symme- 
try of the spacetime. We now specify it further and as- 
sume that it is close to a flat Friedmann-Lemaitre space- 
time (FL), which we consider as our background space- 
time. The homogeneity and isotropy assumptions imply 
that 



a 2 = -a ab a ab = , tu ab = , V./ = , (26) 
where / is any scalar quantity; in particular 

V„p = V a p = => V a ii = , a a = . (27) 
The background (zero - order) equations are given by [3^ | : 

3H + 3H 2 = 8ttG [V(cf>) - ip 2 ] , (28) 
3H 2 = 8irG [\i\} 2 + V(<j»)] (29) 
i) + ZHip + V'{4>) = & ft + ZH1P 2 = , (30) 

where all variables are a function of cosmic time t only. 



III. GRAVITATIONAL WAVES FROM DENSITY 
PERTURBATIONS: COVARIANT FORMALISM 

A. First order equations 

The study of linear perturbations of a FL background 
are relatively straightforward. Let us begin by defining 
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the first-order gauge-invariant (FOGI) variables corre- 
sponding respectively to the spatial fluctuations in the 
energy density, expansion rate and spatial curvature: 

X a = V„/j , 

Z a = v a e, 

C a = a 3 V a R. (31) 

The quantities are FOGI because they vanish exactly in 
the background FL spacetime [38j, [40( . It turns out that a 
more suitable quantity for describing density fluctuations 
is the co-moving gradient of the energy density: 

V a = ~X a , (32) 
M 

where the ratio X a /fi allows one to evaluate the mag- 
nitude of the energy density perturbation relative to its 
background value and the scale factor a guarantees that 
it is dimensionless and co-moving. 

These quantities exactly characterize the inhomogene- 
ity of any fluid; however we specifically want to charac- 
terize the inhomogeneity of the scalar field: this cannot 
be done using the spatial gradient V Q </> because it iden- 
tically vanishes in any space-time by virtue of our choice 
of 4- velocity field u a . It follows that in our approach the 
inhomogeneities in the matter field are completely incor- 
porated in the spatial variation of the momentum den- 
sity: V„^, so it makes sense to define the dimensionless 
gradient 

* a = ^V a ?/> , (33) 

which is related to T> a by 

ib 2 

V a = "—^ a = 7 * a , (34) 
A* 

where we have used Eq. (fT8|) and 7 is given by Eq. (|23|) : 
comparing Eq. (j3U|) and Eq. (fP3"j) we see that ^f a is propor- 
tional to the acceleration: it is a gauge-invariant measure 
of the spatial variation of proper time along the flow lines 
of u a between two surfaces 4> =const. (see Ref. 33]). The 
set of linearized equations satisfied by the FOGI variables 
consists of the evolution equations 



X a = -AHX a ~i, 2 Z a , (35) 

Z a = -3HZ a - AnGX a + V a div ii, (36) 

Vab - V( a u b ) = -2Ha ab - E abl (37) 

E ab - curl H ab = -4irGip 2 a ab - QE abl (38) 

H ab + cur\E ab = -3HH ab ; (39) 



and the constraints 



= ^fX a -dwE a , (40) 

= §Z Q -diva a , (41) 

= div# a , (42) 

= Hob - curl a ab , (43) 

= curlX a , (44) 

= curl Z a . (45) 

The curl operator is defined by cur\ip ab = (curl^) Q f, = 



£ c d(a^ c ' l Pb} wnere e abc is the completely antisymmetric 
tensor with respect to the spatial section defined by 
(bed = e a bcdu a , tabed being the volume antisymmetric 
tensor such that e i23 = \J— g- The divergence div 
of a rank n tensor is a rank n — 1 tensor defined by 

(div VOu. ..i„_i = V l ™W..i„- 

Because the background is homogeneous and isotropic, 
each FOGI vector may be uniquely split into a curl-free 
and divergence-free part, usually referred to as scalar and 
vector parts respectively, which we write as 

V a = Vs a + V Va , (46) 

where curlK a = and divK = 0. Similarly, any tensor 
may be invariantly split into scalar, vector and tensor 
parts: 

Tab = Ts ab + T Vab + T Tab (47) 

where curlT S(lb = , div divT v = and divT T(1 = 0. It 
follows therefore that in the above equations we can sep- 
arately equate scalar, vector and tensor parts and obtain 
equations that independently characterize the evolution 
of each type of perturbation. In the case of a scalar field, 
the vorticity is exactly zero, so there is no vector contri- 
bution to the perturbations. 

Let us now concentrate on scalar perturbations at lin- 
ear order. It is clear from the above discussion that pure 
scalar modes are characterized by the vanishing of the 
magnetic part of the Weyl tensor: H ab — 0, so the above 
set of equations reduce to a set of two coupled differential 



equations for X a and Z a : 

X a + 4HX a = -i> 2 Z a (48) 

Z a + 3HZ a = -4nGX a -yj- 2 V 2 X a , (49) 

and a set of coupled evolution and constraint equations 
that determine the other variables 

dab = -ip- 2 V {a X b) - 2Ha ab - E ab . (50) 

E ab = -AttG i> 2 a ab - 3HE ab , (51) 

di V( j a = |Z a , (52) 

curl a ab = 0, (53) 

d\vE a = ^fX a . (54) 
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B. Gravitational waves from density perturbations 

The preceding discussion deals with first-order vari- 
ables and their behavior at linear order. It is important 
to keep in mind that we were able to set H ab = only 
because pure scalar perturbations in the absence of vor- 
ticity implies that curl<7 a o = at first order. The van- 
ishing of the magnetic part then follows from equation 
P3"| . However, at second order curl a ab 7^ 0. We denote 
the non-vanishing contribution at second order by [2~0 | 

S ao = curl a ab . 

The new variable is second-order and g aug e-invariant 
(SOGI), as it vanishes at all lower orders [33|. It should 
be noted that the new variable is just the magnetic part 
of the Weyl tensor subject to the conditions mentioned 
above i.e. 

S a b = H a b\ u =o. (55) 

We are interested in the properties inherited by the new 
variable from the magnetic part of the Weyl tensor. In 
particular, it can be shown that the new variable is trans- 
verse and traceless at this order and is thus a description 
of gravitational waves. It should be stressed that in full 
generality, there are tensorial modes even at first order. 
By assuming that there are none, we explore a particular 
subset in the space of solutions. From the "iterative res- 
olution" point of view, this means that we constrain the 
equations in order to focus on second order GWs sourced 
by terms quadratic in scalar perturbations. In doing so, 
we artificially switch off GW perturbations at first order. 

C. Propagation equation 

The propagation of the new second-order variable now 
needs to be investigated using a covariant set of equa- 
tions that are linearized to second-order about FL. We 
make use of Eqs. (|20| . (|2Tj) and the following evolution 
equations which are up to second order in magnitude; 

E ab = -QE ab + cur\Z ab -4irG^ 2 (j ab + 3<T c{a E b) {^) 
t ab = -eS ao - curl E ab - 2e cd{a u c E b) d , (57) 

together with the constraint 

^ = -^V^ = -^div^. (58) 

Unlike at first-order, where the splitting of tensors into 
their scalar, vector and tensor parts is possible, at second 
order this can only be achieved for SOGI variables. 

We may isolate the tensorial part of the equations by 
decoupling S a fc: since it is divergence free it is already a 
pure tensor mode, whereas E ab is not. The wave equa- 
tion for the gravitational wave contribution can be found 
by first taking the time derivative of (|57p and making ap- 
propriate substitutions using the evolution equations and 



keeping terms up to second order. The wave equation for 
S a b then reads: 

Sa -V 2 Sa + 7i/S ao + (12i/ 2 -167rG^ 2 )Sa 6 - Sab (59) 

where the source is given by the cross-product of the 
electric- Weyl curvature and its divergence (or accelera- 
tion): 

S ab = - [2u e V e + 16H - l5Hc 2 ] 

(w ecd<aVdiV ^)' (60) 

To obtain this, we have used the fact that with a flat 
background space-time 

curlcurlTafc = -V 2 T ah + ^V (a div T b) 

and used the commutation relation 

(curl T ab )- = curl f ab - curl T ab ~ e cd{a a ec D^T b) d 

+e cd{a [u c f b d + ^Qu c T b) d }. 

We have also used Eqs. (pM)) and (f2"S"l) to eliminate ip/i/) 
from the source term. It can also be shown that S ab is 
transverse, illustrating that Eq. (|59)) represents the grav- 
itational wave contribution at second order. Note that 
this is a local description of gravitational waves, in con- 
trast to the non-local extraction of tensor modes by pro- 
jection in Fourier space. Since S Q0 contains exactly the 
correct number of degrees of freedom possible in GW, 
any other variable we may choose to describe GW must 
be related by quadrature, making this a suitable master 
variable. The situation is analogous to the description 
of electromagnetic waves: Should we use the vector po- 
tential, the electric field, or the magnetic field for their 
description? Mathematically it doesn't matter of course 
- each variable obeys a wave equation and the others 
are related by quadrature. Physically, however, it's the 
electric and magnetic fields which drive charged particles 
through the Lorentz force equation - the electromagnetic 
analogue of the geodesic deviation equation. 

In order to express the gravitational wave equation in 
Fourier space, we define our normalised tensor harmonics 
as 

cab 



where £ ab is the polarization tensor, which satisfies 
the (background) tensor Helmholtz equation: \7~Q ab = 
— (q 2 /a 2 )Q ab . As q a is required to satisfy q a u a = in the 
background, it can thus be identified with a 3- vector and 
will subsequently written in bold when necessary. We 
denote harmonics of the opposite polarization with an 
overbar. Amplitudes of S ao may be extracted via 

E(M) = [ d 3 k[Z ab ( Xl t)Q* ab {k lX )], (62) 
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with an analogous formula for the opposite parity. This 
implies that our original variable may be reconstructed 
from 



. (63) 

The same relations hold for any transverse tensor. Hence, 
our wave equation in Fourier space is 

S"(fc, ij) + 6WE'(fe, rf) + [k 2 + 12H 2 - 16ttG^ 2 ]£(A;, rj) 

= S(k, V ), (64) 

with an identical equation for the opposite polarization. 
We have converted to conformal time rj, where a prime 
denotes derivatives with respect to rj, and we have de- 
fined the conformal Hubble parameter as TL = a' J a. The 
source term is composed of a cross-product of the electric 
part of the Weyl tensor and its divergence. At first-order, 
the electric Weyl tensor is a pure scalar mode, and can 
therefore be expanded in terms of scalar harmonics. To 
define these, let = e*** / (2n) 3 / 2 , be a solution to 
the Helmholtz equation: V 2 Q^ = -{q 2 / a 2 )Q {s) . Begin- 



ning with this basis, it is possible to derive vectorial and 
(PSTF) tcnsorial harmonics by taking successive spatial 
derivatives as follows: 



S afc = J d 3 fc [E(fc, t)Q ab (k, x) + E(fc, t)Q ab (k, x)} . Qi s) 



)«:» 

' at} 



a 



Qa % 



Q 



(65) 



(66) 



This symmetric tensor has the additional property 
Q a q b Qab = ~(2g ,4 /3a 2 )Q^ s ^. Using this representation we 
can express our source in Eq. (|64|) in terms of a convo- 
lution in Fourier space, by expanding the electric Weyl 
tensor as 



E(q,n) 



(67) 



Then, the right hand side of Eq. (|60p expressed in confor- 
mal time, accompanied by appropriate Fourier decompo- 
sition of each term and making use of the normalization 
condition for orthonormal basis, yields: 



S(k,r,) = J d 3 qA(q, k){2[E(q, r,)E(k ~ q, V )}' + (16 - 15c 2 )H E(q, V )E(k - q,r,)} 



where 



A(q,k) = 



67rGa 3 ip 2 



cd{ aq b) q d (k c -q c )\k-q\ 2 C b (k), 



(68) 



(69) 



with a similar expression for the other polarization. 

In principle we can now solve for the gravitational wave contribution E b, and calculate the power spectrum of 
gravitational waves today. For this however, we need initial conditions for the electric Weyl tensor (or, alternatively 

IV. GRAVITATIONAL WAVES FROM DENSITY PERTURBATIONS: COORDINATE BASED 

APPROACH 



In this formalism, we consider perturbations around a FL universe with Euclidean spatial sections and expand the 
metric as 



ds 2 = a 2 (r)) {-(1 + 2A)dr 1 2 + 2u i dx i dr] + [(1 + 2C)<5 y - + h tJ ] dx l dx j ) 

I 



(70) 



where 77 is the conformal time and a the scale factor. We 
perform a scalar-vector-tensor decomposition as 



uo t = DiB + Bi 



and 



h — 2f 



where Bi, £j are transverse (E>i£- % 
traceless and transverse (£\ — DiB 1 , 



(71) 



(72) 



2DiDj£ , 

= DiB 1 ), and £y is 
= 0). Latin indices 



i,j,k... are lowered by use of the spatial metric, e.g. B l — 
^Bj. We fix the gauge and work in the Newtonian 
gauge defined by Bi = £ = B = so that $ = A and VP = 
— C are the two Bardeen potentials. As in the previous 
sections, we assume that the matter content is a scalar 
field 4> that can be split into background and perturbation 
contributions: (f> = (j>(rj) + 5<t>(r),x). The gauge invariant 



7 



scalar field perturbation can be defined by 



Q = S<f»-<f/ 



C 

n ' 



where TL = a' /a = aH. We denote the field perturbation 
in Newtonian gauge by x so that Q = \ + (cf)' /H)^ . 
Introducing 



Once the initial conditions are set, solving Eq. (|79|) will 
give the evolution of v k {rf) during inflation, from which 
(73) ®k{v) ancl T^-kiv) can be deduced, using the previous ex- 
pressions. 

Defining the power spectrum as 

(TZ k n* kl ) = ^-Vn(k)S^(k - k'), (82) 



£ ~ 2 n 



(74) one easily finds that 



the equation of state ([23)1 takes the form 7 — w + 1 = 
2e/3. We thus have two expansions: one concerning the 
perturbation of the metric and the other in the slow-roll 
parameter e. 



Scalar modes 



= 2^ 



(83) 



Note also that z and e are related by the simple relation 
V^rGz = ay/e, (84) 

so that 



Focusing on scalar modes at first order in the pertur- 
bation, it is convenient to introduce 



and 



v = aQ 



H 



in terms of which the action ([T]) takes the form 
1 



Sscal = r / d Xdr ) 



(v'Y - (d iV y + 



(75) 



(76) 



(77) 



when expanded to second order in the perturbations. It is 
the action of a canonical scalar field with effective square 
mass m 2 = —z"/z. v is the canonical variable that must 
be quantized [4l|. It is decomposed as follows 

d 3 fc 



V(X,T]) 



(2tt) 3 / 2 



v k (r])e lk - x a k + v* k (r])t 



ife.sc't 



Here v k is solution of the Klein-Gordon equation 



k 2 )v k = 



(78) 



(79) 



and the annihilation and creation operators satisfy the 
commutation relation, [afc,a^./] = S(k — k'). We define 
the free vacuum state by the requirement fife|0) =0 for 
all k. 

From the Einstein equation, one can get the expression 
for the Bardeen potential (recalling that = $) 



A$ = 4irG- 



a 2 $ 



= AttGzv 



(80) 



n \z) 1 \n 

and for the curvature perturbation in comoving gauge 
K = -v/z. (81) 



X =Q--$=Q- 

a 



4ttG 



4>. 



(85) 



B. Gravitational waves at linear order 

At first order, the tensor modes are gauge invariant 
and their propagation equation is given by 



+ mi' l3 - ASij = 



(86) 



since a minimally coupled scalar field has no anisotropic 
stress. Defining the reduced variable 



(87) 



the action (Q]) takes the form 



Stens = \ J ^X&'n 



when expanded tosecond order. Developing £y , and sim- 
ilarly Hij, in Fourier space: 



A=+,x 



d 3 fc 



£xeUk)e 



ik.x 



( 27r )3/2 a „ 



(89) 



where £y is the polarization tensor, the action ([88]) takes 
the form of the action for two canonical scalar fields with 
effective square mass m 2 = —a" /a 



(Ma) - (dkVxY + -Hi 
a 



(90) 

If one considers the basis (ei,e2) of the 2 dimensional 
space orthogonal to k then = (e}e^ — e 2 e 2 )S^_ + (eje 2 + 
e 2 e))Sl 
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/<a arc the two degrees of freedom that must be quan- 
tized [4l| and we expand them as 



d 3 k 



(2tt) 3 / 2 



fc.A 



Hk is solution of the Klein-Gordon equation 



/4' + ^ - — U* = o 



4(fe)(91) 



(92) 



where we have dropped the polarization subscript. The 
annihilation and creation operators satisfy the com- 
mutation relations, [frfc,A>frfc< v] = 8(k — k')5\\i and 
[Sfej^L a! = 0- We define the free vacuum state by the 
requirement 6fc, a 1 0} = for all fc and A. 
Defining the power spectrum as 



(£k,\Et>,x>) = ^V T (k)5 {3) (k - k')6x X , 



one easily finds that 



V T (k) = 16ttG o ., 



fc 3 //fc 2 



(93) 



(94) 



where the two polarizations have the same contribution. 



C. Gravitational waves from density perturbations 



At second order, we split the tensor perturbation as 
£ij = By + £jj /2. The evolution equations of £-^ is 



similar to Eq. (|86[) . but inherits a source term quadratic 
in the first order perturbation variables and from the 
transverse tracefree (TT) part of the stress-energy tensor 



a 2 [t;] TT = 7 ^ [d jX d P xf T - 
It follows that the propagation equation is 

£S ] " + *HgW - A£f = 



(95) 



(96) 



where Sjj T is a TT tensor that is quadratic in the first 
order perturbation variables. 



Working in Fourier space, the TT part of any tensor 
can easily be extracted by means of the projection oper- 
ator 



-L^j (fc) — 6%^ 



(97) 



where fc = k % jk (note that _Ly (fc) is not analytic in k 
and is a non-local operator) from which we get 



S??(k, V ) 



a I b 
-i -Lj 



ab 



S a b(k,r)) 



EE P^(k)S ab (k, V ). 



(98) 



The source term is now obtained as the TT-projection 
of the second order Einstein tensor quadratic in the first 
order variables and of the stress-energy tensor 



S a b — S 



(2) 



I o( 2 ) I q{ 
ab,SS t ° a b,ST ' °ab,TT 



j(2) 



(99) 



The three terms respectively indicate terms involving 
products of first order scalar quantities, first order scalar 
and tensor quantities and first order tensor quantities. 
The explicit form of the first term is 



s: 



TT 



TT 



(100) 



The first term was considered in Ref. [43[ and the second 
term was shown to be the dominant contribution for the 
production of gravitational waves during preheating [12] • 
In Fourier space, it is given by 



S* (2) - -4 



J d 3 qqbqaHq,v)®(k- q,v) 
4ttG J d 3 qq b q a x(Q,il)x(^- Q,V) 



(.101) 



(2) _____ 

fj,\-' (x,r)) can be decomposed as in Eq. using the 



same definition (]8T[) at any order. The two polarizations 
evolve according to 



(2)1/ 

Ma 



fc 2 - — 



/'a 



(2) 



2 a 



pafc c( 2 ) Jj 
r ij °ab,SS t \ ■ 



(102) 



Since the polarization tensor is a TT tensor, it is obvious 



that I']''' [{ - - ab 



so that 



~y= s x J d 3 9ft<Zj ;[${q,ri)$(q-k,r)) + 47rGx(q,vMq-k,r))]- (103) 



From the equation (|102p , we deduce that the source term derives from an interaction Lagrangian of the form 

Stat = / d v d 3 x-^L= [d^dj® + AwGdixdjx] //■ 
J v ottG 



(104) 
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It describes a two-scalars graviton interaction. In full 
generality the interaction term would also include, at 
lowest order, cubic terms of three scalars, two gravitons- 
scalar and three gravitons. They respectively correspond 
to second order scalar-scalar modes generated from grav- 
itational waves and second order tensor modes. As em- 
phasized previously, we do not consider these interactions 
here. 



V. COMPARISON OF THE TWO FORMALISMS 

Before going further it is instructive to compare the 
two formalisms and understand how they relate to each 
other. Note that we go beyond Ref. [35[ , where a compar- 
ison of the variables was made at linear order. Here we 
investigate how the equations map to each other and ex- 
tend the discussion to second order for the tensor sector. 
At the background level the scale factors a and expan- 
sion rates H introduced in each formalism agree, which 
explains why we made use of the same notation. 

The perturbations of the metric around FL space-time 
has been split into a first-order and a second-order part 
according to 



(105) 



We make a similar decomposition for the quantities used 
in the 1 + 3 covariant formalism. As long as we are in- 
terested in the gravitational wave sector, we only need to 
consider the four-velocity of the perfect fluid describing 
the matter content of the universe which we decompose 
as 



a 



Its spatial components are decomposed as 

V 1 = &v + v l , 



(106) 



(107) 



V 1 being the vector degree of freedom and V the scalar 
degree of freedom. As V^ 1 has only three independent de- 
grees of freedom since satisfies u^u^ = —1, its tempo- 
ral component is linked to other perturbation variables. 
We assume that the fluid has no vorticity (V 1 =0), as it 
is the case for the scalar fluid we have in mind and con- 
sequently we will also drop the vectorial perturbations 
(Si = 0). 



A. Matching at linear order 

At first order, the spatial components of the shear, 
acceleration and expansion are respectively given by 

a^=a(d {i d j} V^+£^') , (108) 
u^ = di^+HV^+V^') , (109) 



(110) 



The electric and magnetic part of the Weyl tensor take 
the form 

1 



e\;' =d {i d j) ^-^{e': j + A£ i3j 



H 



(i) 

y 



(111) 



(112) 



Note that r}ku is the completely antisymmetric tensor 
normalized such that 77123 — 1, which differs from E a bc- 
We deduce from the last expression that 



curl 



1 

2a 



(curl £«" 



curl AS (1) 



A 13 ) 

where we have used simpler notation by writing (curl S)ij 
as curlfij. We also note that the derivative along of 
a tensor T of rank (n, m), vanishing in the background, 
takes the form 

fil h. in 3m = d * T XX +( U - m ) HT \i n i m ( 114 ) 

at first order, or alternatively 



Ox — 



= d t r 



(115) 



Again, recall that a dot refers to a derivative along u^. 
Indeed at first order, it reduces to a derivative with re- 
spect to the cosmic time but this does not generalize to 
second order. 

Now, Eq. (f39|) can be recast a 



a- 2 (a^Hij)' + curl Ey + H H tJ = . (116) 

Using the expressions (|HHI112j) for the geometric quan- 
tities, this equation takes the form 



curl 



AS,, 



= 0. 



(117) 



Similarly Eq. ([59")) can be recast as 



(a 2 H ab 



3H 



.a 2 H a b 



so that it reduces at first order to 



2{H 2 + H)H ab -W 2 H ab = 0, 
(118) 



curl 



1 



2a 2 ^ 



-^r ( + 2HS 



(i)' 

ij 



AS 



(i) 



= 0. 



(119) 



Thus, Eq. (TIT9"|) maps to Eq. ([SB) with the identifi- 
cation pi2[) , if there is no vector modes. This can be 
understood from the fact that in the Bardeen formal- 
ism, Eq. (|86[) is obtained from the Einstein equation as 

curl [curl G,j] = 0. 

In the case where there are vector modes, Eq. (| 1 1 2[) 
has to be replaced by 



H 



(i) _ 



= {cuv\S^%+ l - m(l d k d 0) S^ n 



and H a i, is no longer a description of the GW, i.e. directly 
related to the TT part of the spacetime metric and the 
matching is not valid anymore. 
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B. Matching at second order 

At second order, the matching is much more intricate 
mainly because the derivative along w M does not match 
with the derivative respect to cosmic time any more. 

Let us introduce the short hand notation 

(X x 1%. = Vki(iX k Y j} 1 (120) 



for any tensors X k and Y lm . If Y lm = d l d m Z, or X k = 
d k W, we also use the short-hand notation Y = ddZ X = 
dW. 

Among the terms quadratic in first-order perturba- 
tions, those involving a first-order tensorial perturbation 
can be omitted, as we are only interested in second-order 
effects sourced by scalar contributions. At second order, 
the geometric quantities of interest read 



H 



(2) 



(curl£ (2) 



(curl - 4 (dV {1) x dd^> {1 

(curl£( 2 >") + (curl AS {2 



1 

2^ 



(<9$ (1) x ,9<9<I> (1) ) +H (dV w x <9<9$ (i; 
- (dV {1) x <9«9$ (1) ' N 



(121) 



(122) 



From the latter expression, we remark that H\a has a 
term quadratic in first-order perturbations involving 
and This terms arise from a difference between 

the two formalisms related to the fact that geometric 
quantities, such as Hij Eij etc., live on the physical space- 
time, whereas in perturbation theory, any perturbation 
variable at any order, such as etc., are fields 

propagating on the background space-time. 

It follows that the splitting into tensor, vector and 
scalar modes is different. In the covariant formalism, 
the splitting refers to the fluid on the physical space- 
time, whereas in perturbation theory it refers to the co- 
moving fluid of the background solution. Indeed, this 
difference only shows up at second order as the magnetic 
Weyl tensor vanishes in the background. The one to one 
correspondence at first order between equations of both 
formalisms disappears, as the second order equations of 



the covariant formalism contain the dynamics of the first 
order quantities. 

When keeping terms contributing to the second order, 
Eq. ([55)) has an additionnal source term and reads 

H ab + curl E ab + 3HH ab = -2e cd(a u c E b f (123) 

If first order tensorial perturbations are neglected then 
H ab vanishes at first order and Eq. (|114p still holds when 
applied to H ab . Thus Eq. (|123p can be recast as 

(a 2 H ab \ U fi 

± P- + curl Eat + ~H ab = -2e cd(a u c EJ . (124) 

Substituting the geometric quantities for their expres- 
sions at second order, and making use of Eq. (|115[1 to 
handle the derivatives, Eq. (|116[) reads at second order 



curl 



1 (p{2)ll 

2a r« 



2HB[f - A£g ] 



(a$ (1) x - (dv {1) x - n (ay (1) x 



Using the momentum and constraint equation (|41|) at first order 



(125) 



(126) 



and the background equation H! — 7i 2 = — 47rG/i(l + w)a 2 , that we deduce from the Raychaudhuri equation and the 
Gauss-Codacci equation at first order, we can link it to Eq. as it then reads 



— curl 

a 



lie 



5(2)" 



2H£[f - AEg ] 



= —curl 

a 



2dM 1) d J <& (1) + SnGa 2 {ix + P)d i V^d j V w 



(127) 



When applied to a scalar field, this is exactly the gravitational wave propagation equation 
term (TTIJOl . 

I 



with the source 



Discussion perturbations. This extends the work of Ref. [35j which 



In conclusion, we have matched both the perturbation 
variables and equations at first and second order in the 



11 



considered the linear case, and has not been previously 
investigated. 

Even though we restrict to the tensor sector, this com- 
parison is instructive and illustrates the difference of ap- 
proach between the two formalisms, in a clearer way than 
at first order. In the Bardeen approach, all perturbation 
variables live on the unperturbed spacetime. At each or- 
der, we write exact equations for an approximate space- 
time. In particular, this implies that the time derivatives 
are derivative with respect to the cosmic time of the back- 
ground spacetime. In the covariant approach, one derives 
an exact set of equations (assuming no perturbation to 
start with). These exact equations are then solved it- 
eratively starting from a background solution which as- 
sumes some symmetries. The time derivative is defined 
in terms of the flow vector as u a V a . Indeed, at first or- 
der for scalars, this derivative matches exactly with the 
derivative with respect to the background cosmic time. 
At second order, this is no longer the case. First the 
flow vector at first order does not coincide with its back- 
ground value. This implies a (first-order) difference be- 
tween the two time derivatives which must be taken into 
account. Then, the geometric quantities, such as E%j 
etc., "live" on the physical space-time, whereas in per- 
turbation theory, any perturbation variable at any order, 
such as etc., live on the background space- 

(2) 

time. This explains why e.g. has a term quadratic 

in first-order perturbations involving V' 1 ) and 

The master variables and corresponding wave equa- 
tions in both formalisms are also different in nature. In 
the metric approach the wave equation with source is 
defined non-locally in Fourier space; in the covariant ap- 
proach, we are able to derive a local tensorial wave equa- 
tion which, because it is divergence-free, represents the 
gravitational wave contribution. Of course, we can make 
a non-local decomposition in Fourier space as required. 
Furthermore, on one hand the TT part of the metric in 
a particular gauge is a perturbative approach used to de- 
scribe GW, and this tells us the shear of spatial lengths 
with respect to a homogenous and isotropic background, 
referring implicitly to a hypothetical set of averaged ob- 
servers. On the other hand, the covariant description 
using H a b which is built out of the Weyl tensor and the 
comoving observer's velocity, directly describes the dy- 
namically free part of the gravitational field [39[ (up to 
second-order when rotation is zero) as seen by the true 
comoving observers. This is part of the dynamic space- 
time curvature which directly induces the motion of test 
particles through the geodesic deviation equation, and it 
accounts for effects due to the non-homogenous comoving 
fluid velocity. 

There is one more difference between the two for- 
malisms, concerning the initial conditions. In the 
Bardeen approach, as we recalled in section HVi there 
is a natural way to set up the initial conditions on sub- 
Hubble scales by identifying canonical variables, both for 
the scalar and tensor modes, and promoting them to the 



status of quantum operators. In the covariant formalism 
such variables have not been constructed in full general- 
ity (see however Ref. [44| for a proposal). Consequently 
this sets limitations to this formalism since it cannot ac- 
count for both the evolution and the initial conditions at 
the same time. 



VI. ILLUSTRATION: SLOW-ROLL INFLATION 

A. Slow-roll inflation 

In this section, we focus on the case of a single slow- 
rolling scalar field and we introduce the slow-roll param- 
eters 



3^ 2 
27' 



Htp 



TT7- (128) 



Using the Friedmann equations (|28M29[) . these parameters 
can be expressed in terms of the Hubble parameter as 



1 



H'(4>) 



4ttG 
Interestingly Eq 

which implies 



, 1 H"{4>) 

H{4>) J ' AttG H(4>) ' 

takes the form 



H 2 



1 



a 
a 



(129) 



(130) 



(131) 



The equation of state and the sound speed of the equiv- 
alent scalar field are thus given by 



w 



-1 



(132) 



The evolution equations for e and d show that e and 5 
are of order 2 in the slow-roll parameters so that at first 
order in the slow-roll parameters, they can be considered 
constant. Using the definition of the conformal time and 
integrating it by parts, one gets 



a(rj) 



1 1 



Hi] 1 — s 

assuming e is constant, from which it follows that 
H = aH = — (l + e) + 0(2), 



(133) 



(134) 



where r\ varies between — oo and 0. This implies that 



a 
a 



2 + 3e 



2 + 6e - 36 



The general solution of Eq. ({T9l is 



(135) 



Vk 



= V-7Tf 1 /A aHjPi-kr]) +c 2 H& (-k v ) , (136) 
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with |ci| 2 — | C2 1 2 = 1, where Hjp and are Hankel 
functions of first and second kind and z^ = 3/2 + 2e — S. 
Among this family of solutions, it is natural to choose the 
one with c 2 = which contains only positive frequen- 
cies [4l|. It follows that the solution with these initial 
conditions is 



Similarly as for the scalar mode, we obtain 



V T (k) 



16 H 2 

~W 2 l 



-2s 



(141) 



On super-Hubble scales, \krj\ <C 1, we have 

v k -> 2^ 3 / 2 r(^)/r(3/2)(2fc)- 1 / 2 (-fc7 ? )- l/+1 / 2 . 



(137) 



Now, using Eq. (|133|) to express rj and Eq. ([84]) to replace 
z in expression (|83p . we find that 



Vn(k) 



1 H 2 



7T M 2 e 



2^-3/2 



r(3/2) 

k 

aH 



2 / , N -2i/+l 

2/ — — 



-2^+3 



(138) 



where we have set M 2 



G 1 . At lowest order in the 
slow-roll parameter, it reduces to 



^ 1 H 2 ( k \ 



25-is 



(139) 



The evolution of the gravitational waves at linear order 
are dictated by the same equation but with vr = 3/2 + e, 
so that 



(140) 



B. Gravitational waves at second order 

The couplings between scalar and tensor modes at sec- 
ond order imply that the second order variables can be 
expanded as 



1 



,(2) 



n 



(2) 

as 



and a similar expansion for £, where, e.g., TZ^£ stands for 
the second order scalar modes induced by the coupling 
of first order scalar and tensor modes etc. The deviation 
from Gaussianity at the time rj of the end of inflation can 
be characterized by a series of coefficients defined 
for example as 

^(M) - J 5 3 (k 1 + k 2 -k)£(k u rj)£{k 2 ,rj) 



f£i ££ {k,k 1 ,k 2 ,r))d 3 k 1 d 3 k 2 . (142) 

These six coefficients appear in different combinations in 
the connected part of the 3-point correlation function of 
1Z and £. For instance 



(£ kl K k2 n k3 ) c = 2f^ K (k 1 ,k 2 ,k 3 )P n (k 3 )P n (k 2 ) 



f&? 1i (k 1 ,k2,k 3 )P R .(k 3 )Pe(k 1 ) a a (fci + fc2 + fc 3 ) (N:5) 
I 



and /j 



NL 



is the standard /nl parameter. One in the slow-roll parameter to 



can easily check that ( r R-k 1 T^k 2 '^-k3)c involves 
/nl RR . (ffc 1 ^fc 2 ^fc 3 ) c involves / N £ n and f } 
and (£k 1 £k 2 £k 3 )c involves f } 



£,££ 
NL 



K,££ 
'NL 



C. Expression for /^l^ 



From our analysis, we can give the expression of /nl^ • 
Starting from the fact that -ell = $(1 + e) + & /H and 
from the expression (|85|) . we get that $ ~ — elZ — $'/7Y 

or $ = -erjf Zrdrj, and VI^G X - & /H]. It 

follows that the source term <| 100|) reduces at lowest order 



£TT = 4 [ ed . nd . n 



TT 



The interaction Lagrangian is thus given by 



= / d?yd 3 a; 



4a 



ediRSjUn* 3 , 



which reduces to 

S int = ( dr)d 3 x2d i vd j v£ ij . 



(144) 



(145) 



This is the same expression as obtained in Ref. @. 

In full generality, during inflation, we should use the 
"in-in" formalism to compute any correlation function of 
the interacting fields. As was shown explicitly in Ref. [45[ 
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for a self- interacting field and more generally in Ref. [47[ , 
the quantum computation agrees with the classical one 
on super-Hubble scales at lowest order. Note however 
that both computations may differ (see Ref. [26| versus 
Ref. [9() due to the fact that in the classical approach 
the change in vacuum is ignored. The difference does 
not affect the order of magnitude but the geometric k- 
dependence. In order to get an order of magnitude, we 
thus restrict our analysis here to the classical description. 
This description is also valid when considering the post- 
inflationary era. 

In the classical approach, we can solve Eq. p03[) 
by mean of a Green function. Since the two in- 
dependent solutions of the homogeneous equation are 
V ' —krjH^j/ 2 ^ {—krj), the Wronskian of which is Ai/{-nk), 
the Green function is given by 



It follows that the expression of the second order tensor 
perturbation is given by 



" dr,'4=eG(k,r,,r,') 



(2^)3/2 J^-' 

d 3 q (Wy) K,(»/)ft fc -,(»/)<147) 



G(k, T],T]') 



Hg>(-kT,)H&>(- 



■krf) 



H^(-kr/)Hj^{-krj)\. (146) We thus obtain 



-l £ 



a(v) J- 



drfatf) G{k, 77, rf) {qudij^ (*)) K qi (r/)^ 2 (rj') . (148) 



r 



If we want to estimate Eq. (|143[) in the squeezed 
limit k\ <C k2, k^ the contribution coming from the 
term involving f N x ^ ln (k,qi,q2,r)) can be computed by 



integration times. From the Wick theorem, this corre- 
lator reduces to n(q,7]')TZ*(q,r]")n(\k - q\,rf)n*(\k - 
q\,rj")S(k — k')[8(q — p) + S(k — q — p)] and because 



use of the super-Hubble limit of the Green function fc l £y = the two terms give the same geometric factor. 



\G(k,r,,rf)\ 



2v T 











(C 











This contribu- 



(feifey A ) * 3 (*1 



tion will be proportional to irp^k 2 

k2 + ks), which is the same order of magnitude as in 
Ref. [9( , but do not have the same geometric dependence 
as it goes like k 2 5 k^ 3 instead. 



D. Orders of magnitude 



Thus, the integration on p is easily done and we can 
factorize 5(k — k'). Now, note that the terms in the 
integral involve only the modulus of q and k — q so that 
it does not depend on the angle ip of q in the plane 
orthogonal to k. This implies that the integration of ip 
will act on a term of cos 2 2<p, sin 2 2ip and cos 2ip sin 2ip 
respectively for ++, x x and +x so that it gives a term 
TrS\y. In conclusion, defining the second order power 
spectrum V^P by 



When we want to estimate xirftEy y(v))i 
we have to evaluate the connected part of 
(^ <J (7 7 ')^ fc _ <J (7 7 ')^;(r ? ")^^_ p (r ? ")), where q and p 
are the two internal momentum and rj and rj" the two 



\{4%£«,y) = ^-v¥\k)6^(k k')6 X y, (149) 



it can be expressed as 



J 



dr/dr/'a(r/)a(r/') e 2 G(fc, V , ij')G* (fc, V , rj 

2 



(27r) 3 7r 2 a 2 

/ d 8 g(? i ^yw(g,»/)^*(g,»/0^(l*-«ll,'/)^*(l*-<ll,'/ / )- 



(150) 
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Setting k ■ q = kq/x, this reduces to 

A- 3 



(27r) 3 7ra 2 



g 6 dg(l- M 2 ) 2 d M 



dr/ a{rf) e G(k, v , rf)H{q, if)H{\h - q\,rf) 



(151) 



after integration over ip which gives a factor 7r (l — /j?) q 4 



We can now take the super-Hubble limit of this expres- 
sion at lowest order in the slow-roll parameters. In order 
to do so, we make use of the super-Hubble limit of the 
Green function given above, and we perform the time 
integral from 1/fc to rj and keep only the leading order 
contribution: 



1 



3 4 2 3 7T 



-it 



(152) 



where, with the defnitions y = q/k and n = k/k, 



F(e,S) 



J (y\n-y\)- 3 -^ +2S y e dy(l-^) 2 df, 



(153) 

is a numerical factor. In this approximation, the ratio 
between the second order power spectrum and the first 
order power spectrum at leading order in the slow-roll 
parameters, is given by: 



i 



2 7 3 4 tt V M, 



(154) 



Indeed there are ultraviolet and infrared divergences 
hidden in F(e,5). We expect the infrared divergence not 
to be relevant for observable quantities due to finite vol- 
ume effects (see for instance Ref. (46|). The ultraviolet 
divergence, on the other hand, has to be carefully dimen- 
sionaly regularized in the context of quantum field theory 
(see e.g. Ref. [47j). 



VII. CONCLUSIONS 

In this article we have investigated the generation of 
gravitational waves due to second order effects during 
inflation. We have considered these effects both in the 
covariant perturbation formalism and in the more stan- 
dard metric based approach. The relation between the 
two formalisms at second-order has been considered and 
we have discussed their relative advantages. This com- 
parison leads to a better understanding of the differences 
in dynamics between the two formalisms. 

As an illustration, we have focused on GW generated 
by the coupling of first order scalar modes. To charac- 
terize this coupling^we have introduced and computed 
the parameter . It enters in the expression of 

(£k 1 T^k 2 T^k 3 )c that was shown to be of order (i?/A/ p ) 4 /e, 
as (7?.fc 1 7?.fc 2 7?.fc3) c . On the other hand the power spec- 
trum of GW remains negligible. 

This shows that the contribution of {£k 1 TZ-k 2 'R-k 3 )c to 
the CMB bispectrum is important to include in order to 
constrain the deviation from Gaussianity, e.g. in order 
to test the consistency relation [48| . 
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